Spin matrix elements in 2D Ising model 
on the finite lattice 



A. I. Bugrij *, O. Lisovyy *' t 



* Bogolyubov Institute for Theoretical Physics 
Metrolohichna str., 14-b, Kyiv-143, 03143, Ukraine 

' Laboratoire de Mathematiques et Physique Theorique CNRS/UMR 6083, 
' Universite de Tours, Pare de Grandmont, 37200 Tours, France 

o ■ 
o 

<N 

(3JT)! We present explicit formulas for all spin matrix elements in the 2D Ising model with the 

nearest neighbor interaction on the finite periodic square lattice. These expressions generalize 
the known results (TJ [3l [6] (coincide with them in the appropriate limits) and fulfill the test 
■ of straightforward transfer matrix calculations for finite N. 
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Abstract 



1 Eigenvalues and eigenvectors of the transfer matrix 

It is well-known (see 0[8]) that the spectrum of 2^ x 2 N transfer matrix, corresponding to Ising 
model on the periodic square lattice, consists of two sets: 



A = {2s) N/2 exp | i (±7 (0) ± j(2tt/N) ±...± 7 (2vr - 2vr/iV))| , (1) 
C^: A = (2 S ) Af / 2 exp|i(± 7 (^/iV)±7(3^/iV)±...±7(2vr-^/iV))|, (2) 



where s = sinh2/C and JC is the Ising coupling constant. The function -y(q) is defined as the 
positive root of the equation 

cosh j(q) = s + s _1 — cos q. 



which is the lattice analog of the relativistic energy dispersion law. The number of minuses in 
(1) is even in ferromagnetic (s > 1) and odd in paramagnetic (0 < s < 1) phase, while the 
^ ■ number of minuses in (2) is even in both phases. The eigenvalues (1) (or (2)) correspond to 

eigenvectors that are odd (resp. even) under spin reflection. 

The notation and terminology, introduced in [6j for the analysis of continuum limit, are also 
very convenient on the lattice. In what follows, odd and even eigenvectors of the Ising transfer 
matrix will be interpreted as multiparticle states from the Ramond and Neveu- Schwartz sector. 
Quasimomenta of R-particles can be equal to j^j (j = 0, 1, . . . , N — 1), while for NS-particles 
they take on the values (j + •>) = 0, 1, . . . , — 1). Each eigenstate consists of particles of 
only one type, and their quasimomenta must be different. 

We will denote by \pi, . . . ,px) ns(R) the normalized eigenstate, containing particles with the 
momenta p\, . . . ,px- Since R-sector in paramagnetic phase contains the state |0)# (one particle 
with zero momentum), it will be convenient to denote NS and R vacua by \®)ns an d The 
goal of the present paper is to find matrix elements ns{pi, ■ ■ ■ ,PK\&\qi, ■ ■ ■ , qL)R of the Ising spin 
a in the described basis of normalized eigenstates. (R-R and NS-NS matrix elements vanish 
due to Z2-symmetry of the model). 
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2 Lattice form factors and scaling limit 



All n-point correlation functions in the Ising model on the cylinder and torus can be easily 
expressed via spin matrix elements. However, known results were obtained in rather inverse 
way. At the first stage, 2-point functions are expressed through the determinants of certain 
Toeplitz matrices with a size that depends on the separation of correlating spins. To extract 
the analytic dependence on the distance from these representations, a lot of further work was 
needed [2] . The final answer O 0] allows to calculate squared form factors on the cylinder (on 
the infinite lattice the above program was realized earlier in [91 [TOl [T2] ) : 
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Here £ = |l — s -4 ] 1 ^ 4 , quasimomenta have discrete R-values and cylindrical parameters v(q) 
are given by 
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sinh (Nj(p)) sinh (N"f(q)) 
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dp sinh 7(5) 
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sin((p + g)/2) 



2tt J cosh •y(q) — cos p 



sin((p-g)/2) 
lncoth(iV7(p)/2) . 



In the thermodynamic limit N — > 00 these parameters vanish (£y — ► 1, 1/(5) 
transforms into the classical formula [10\ [T2] : 
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0) and (3) 
(6) 



where {qj} can take on arbitrary values in the interval [— ir, tt ]. 

In the scaling limit, Ising model on the plane was shown to be equivalent to a relativistic 
quantum field theory with two-particle S'-matrix equal to —1 (see Then it became possible 

to use the results of pQ and to calculate all spin matrix elements: 
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NS(P1, ■ ■ ■ ,PKW\Ql, ■ ■ -,Ql}r 

f({p}\U} 
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Pi- Pi TT H-Qj TT ^>{Pi) + ^{ij 



n Pi ~ Pj TT W ~ jfl TT 

l<i<j<A Kn ' KrJ ' l<i<j<L yH3 > l<i<K 

1<3<L 



Pi ~ % 



(7) 
(8) 



Here uj(q) = \J m 2 + q 2 and the momenta of both type of particles take on arbitrary real values. 
The RHS of (7) is usually written with the factor but it can be removed by a change of 

the basis and will be omitted in what follows. 

Very recently, Fonseca and Zamolodchikov [6] announced and promised to give a proof of a 
similar formula for the scaling limit on the cylinder: 
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The overall factor £y and the function v{q) from the leg factors are determined from the scaling 
limit of (4), (5): 
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2vr 2 J J sinh (/?w(p)) sinh ' 
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%) = - ! 2 dP f q) 2 Incoth^M. (11) 
w vr J p 2 + q 2 + m 2 2 v 7 
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Here, /? denotes the scaled length of the base of the cylinder, NS-momenta are quantized as 
Pj = lj, lj G Z + 77, while for R-momenta we have = ^ and lj G Z. 

3 General formula 

On the level of form factors n s {^W\{q}) R the expression (9) represents nothing new and can even 
be proven rigorously — it is simply a particular case of (3). However, this conjecture gives all 
matrix elements, though only in the scaling limit. Moreover, the structure of this representation 
is so transparent that the lattice generalization immediately suggests itself. To be more precise, 
having taken into account the correspondence between (3) and (9), we propose the following 
general formula for spin matrix elements on finite periodic lattice: 

JL e "(Pi)/ 2 A. P -"te)/2 

NS (pi, . . . , PK W\qi, ■ ■ ■ jQl)r = V^T 11 -==== 11 



■ Pi-Pj ■ qi-qj ■ . 7(Pi)+7(gj) 

x TT Sm 2 TT Sm 2 TT ^ 2 , 12) 

11 . , tCeiHtM 11 ■_ b 7fe)+7(g,) 11 S in^i ' 1 j 

l<i<j<JsT smn 2 l<i<j<L smn 2 l<i<^ 2 

1<J<£ 

where £r, ^(g) are defined by (4) and (5). All previous results can be easily obtained from 
this expression in appropriate limits. However, we have not yet found a rigorous proof of (12). 
Instead, since this formula should hold even on the finite lattice, we have verified it explicitly 
for small N. 

As an illustration, let us consider 3-row Ising chain in the ferromagnetic region of temperature 
parameter (s > 1). In this case NS (R) momenta take on the values ir, n/3, —tt/3 (0, 2ir/3, 
— 27r/3). Each state contains either two particles or no particles at all. Since the integrals (4) 
and (5) can be alternatively written as 

U {R) U (NS) smh 2 U (NS) sinh tW+tM 

>4_ <? P / \ _ , P 



U {R) U (R) smh iMppl U ^s) n (NS) ginh tW+tM > ™ U ( R) ginh tW+tM 
q p q p " p 



then to verify (12) it suffices to prove ten relations: 



2 sinh 22^Wi sinh 7.+72,/3 ginn 2 7,/3+72,/3 



sinh7 27r /3 sinh 7^/3 sinh 2 2t/3 sinh ■ 
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sinh 2£±^Zi sinh T^Vi sinh 2 7,/3+72,/s 



9 sinh 72^/3 sinh 7^/3 sinh 7o+ 2 7r/3 sinh 7,r+ 2 2 ' r/3 
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ws (0 (t|27t/3, -2tt/3) r 2 



12 sinh 72x/3 sinh 7V3 sinh 2^Zi sinh 2i±^Zi s inh 2 tw»+t*W3 ' 

2 _ S i n h 2^Zi sinh 2i±^Z3 

at S ( n/3,n/3\a\®) R - 7o+7 , /3 7o+72ir/3 2 7j/3+72ir/3 , 



12 sinh 72ff/3 sinh 7w/3 sinh sinh ^ sinh" 



ns (0|<t|O, 2V3)„ 2 = jvs (0k|O, -2tt/3> 2 ' 



12 sinh 22!±^ sinh sinh 2i+^ll s inh , 



jvs (~7r/3,7r|cr|0) fl 2 = ns (tt/3, tt|ct|0) = 



12 sinh 22±2^ sinh 22±^Zi s inh j2w/3 s inh 2d^L ' 



4 sinh 7o+7v / 3 sinh 7,r+77r/3 
ats (-tt/3, 7r/3|<r|0, 2vr/3) i? 2 = at S (-tt/3, 7r/3|cr|0, -2ir/3) R 2 = - - — 2 - - 7 q+ „ 



9 sinh 7rr / 3 sinh ■ 

4 sinh 7o+72,r/3 sinh 7 ' ^+72 ' ^ ' ,3 

ivs (-tt/3, 7r|a|27r/3, -2tt/3) r 2 = NS (tt/3, 7r\a\2n/3, -2tt/3) r 2 = 2 , - 7o+ 2 , 

9 sinh 72w / 3 sinh ' Z, ' 

NS (-tt/3, 7t|<7|0, 2tt/3) r 2 = NS {tt/3, tt|ct|0, -2vr/3) fl 2 = ^ , 

Ns(~n/3,ir\<j\0,-2ir/3) R 2 = Ns(ir/3,ir\<?\0,2n/3) R 2 = g ■ 

This indeed can be done — with a little bit cumbersome but straightforward calculation^. We 
have performed a similar check for small N up to N = 4 and we have no doubt in the validity 
of (12) for arbitrary N. The rigorous proof of this formula will complete, in a sense, the study 
of the 2D Ising model in zero field. 
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